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We investigate the effects caused by noncommutativity of the phase space generated by two 
scalar fields that have non-minimal conformal couplings to the background curvature in the FRW 
cosmology. We restrict deformation of the minisuperspace to noncommutativity between the 
scalar fields and between their canonical conjugate momenta. Then, the investigation is carried 
out by means of a comparative analysis of the mathematical properties (supplemented with some 
diagrams) of the time evolution of variables in a classical model and the wave function of the 
universe in a quantum perspective, both in the commutative and noncommutative frames. We 
find that the impose of noncommutativity causes more ability in tuning time solutions of the 
scalar fields and hence, has important implications in the evolution of the universe. We get 
that the noncommutative parameter in the momenta sector is the only responsible parameter 

^^ | for the noncommutative effects in the spatially flat universes. A distinguishing feature of the 

noncommutative solutions of the scalar fields is that they can be simulated with the well-known 
three harmonic oscillators depending on three values of the spatial curvature. Namely the free, 

[~^. ■ forced and damped harmonic oscillators corresponding to the fiat, closed and open universes, 

respectively. In this respect, we call them cosmical oscillators. In particular, in closed universes, 
when the noncommutative parameters are small, the cosmical oscillators have analogous effect with 
the familiar beating effect in the sound phenomena. Some of the special solutions in the classical 
model and the allowed wave functions in the quantum model make bounds on the values of the 
noncommutative parameters. The existence of a non-zero constant potential (i.e. a cosmological 
constant) does not change time evolutions of the scalar fields, but modifies the scale factor. An 
interesting feature of well-behaved solutions of the wave functions is that the functional form of the 
radial part is the same as the commutative ones within a given replacement of constants caused 
by the noncommutative parameters. Furthermore, the Noether theorem has been employed to 
explore the effects of noncommutativity on the underlying symmetries in the commutative frame. 
Two of the six Noether symmetries of spatially flat universes, in general, are retained in the 
noncommutative case, and one out of the three ones in non-flat universes. 
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1 Introduction 

Scalar field theories have become generic playgrounds for building cosmological models related to 
particle physics [1], and more recently, have played very important contributions in the other branches 
of cosmology, e.g., as a powerful tool for providing expanding accelerated universes. Also, they have 
key roles in some models of the early cosmological inflation [2], and are viable favorite candidates 
for dark matter [3]. Indeed, scalar field cosmological models have extensively been studied in the 
literature, see, e.g., Refs. [4] and references therein. 

Unifying theories of interactions, such as the Kaluza-Klein, supergravity and superstring theories, 
usually predicts non-minimal couplings between geometry of space-time and scalar fields. These 
theories mostly are derived from effective actions of string theories [5] or from compactified Kaluza- 
Klein theories in four dimensions [6]. In both cases, resulting models have non-minimal couplings 
between gravity and one (or even more) scalar field. Actually, considering more than one scalar 
field has been viewed as an easier approach to study accelerating models. These cosmological ideas 
have also been employed in models of inflation to describe the early universe [7] or an evolving 
scalar field known as quintessence [8]. Recently, double scalar-tensor theories have been applied in 
extended gravitational theories, and have given successful descriptions of inflationary scenarios [9, 10]. 
Besides, it has been shown [11] that the reduction of a five-dimensional Brans-Dicke gravity into 
four dimensions is equivalent to a 4-metric coupled to two scalar fields, which may account for the 
present accelerated expansion of the universe. Also, effects of additional scalar fields plus independent 
exponential potentials have been considered in the literature [12]. 

In cosmological models, scalar fields usually present degrees of freedom and appear as dynamical 
variables of corresponding minisuperspaces. This point can be viewed as relevance of noncommuta- 
tivity in these models. The proposal of noncommutativity concept between space-time coordinates, 
in the first time, was introduced in 1947 [13]. Thereafter, a mathematical theory, nowadays known as 
the noncommutative geometry (NCG), has begun to take shape based on this concept since 1980 [14]. 
Noncommutativity among space-time coordinates implies noncommutativity among fields as minisu- 
perspace coordinates. Such a kind of noncommutativity has gotten more attention in the literature. 
In the last decade, study and investigation of physical theories in the noncommutative (NC) frames, 
such as the string and M-theory [15, 16], have caused a renewed interest on noncommutativity in the 
classical and quantum perspectives. In particular, a novel interest has been developed in studying 
the NC classical and quantum cosmologies, e.g. Refs. [17] and references therein. Also, deformation 
of the phase space structure has been viewed as an alternative path, in the context of cosmology, 
toward understanding quantum gravity [18]. The influence of noncommutativity has been exam- 
ined by formulation of a version of the NC cosmology in which a deformation of the minisuperspace 
coordinates [19]-[22] or the minisuperspace momenta [23, 24] is required instead of the space-time 
coordinates. From a qualitative point of view, noncommutativity in the minisuperspace coordinates 
(space sector) leads to general effects, however, a non-trivial noncommutativity in the momentum 
sector introduces distinct and instructive effects in the behavior of dynamical variables. 

In this work, we consider homogeneous and isotropic cosmological models based on a particular 
(multi) scalar-tensor gravity theory with two scalar fields that have non-minimal conformal couplings 
to the space-time. The effects of noncommutativity on the phase space, generated by these fields 
plus the scale factor, are investigated. Our approach is to build a NC scenario via a deformation 
conveyed by the Moyal product [15] in a classical and a quantum perspective, where their cosmological 
implications are more appreciated in the classical one. The procedure of quantization is proceeded by 
the Wheeler-DeWitt (WD) equation for a wave function of the universe by the Hamiltonian operator 
of the model via the generalized Dirac quantization. We introduce noncommutativity between the 
scalar fields and between their canonical conjugate momenta, and will pay more attention to the 
outcome of results. Here, we neglect noncommutativity between scalar fields with the scale factor, 
such effects can be found in, e.g., Refs. [22, 24]. Actually we treat the effects of noncommutativity 
by two parameters which are the NC parameters corresponding to the space and momentum sectors. 



Then, we analyze the mathematical properties of solutions and look for relations, including the ranges 
and values, among the NC parameters for which particular or allowed solutions exist. 

The manuscript is organized as follows. In Section 2, we specify our toy model and investigate 
the classical version within the commutative and NC frames. The quantum version of this model, 
by probing the universe wave functions, is considered in Section 3, where the general properties of 
solutions in the commutative and NC frames are compared. In Section 4, we employ the Noether 
theorem and explore the effects of noncommutativity on the underlying symmetries in the commuta- 
tive frame. Conclusions are presented in the last section, and some necessary formulations have been 
furnished in the appendix. 

2 The Classical Model 

We consider a classical model consisting of a cosmological system presented by a four dimensional 
action with two non-minimally scalar fields coupled to gravity in a Friedman-Robertson-Walker 
(FRW) universe. To specify the NC effects of the model, we first treat the commutative version and 
then consider the NC one in the following subsections. 

2.1 The Commutative Phase Space 

A general action for non-minimally coupled double scalar fields can be written as [9] 

A = J ' v 7 ^ \f{^ i,)R - l -g^ [A(<l>, VOV^V^ + B(4>, V>)V^V„Y>] ~ V{</>, V>)} d 4 x , (1) 

where g is the determinant of the metric g^ u , R is the Ricci scalar, F((f>, ip) and V((f>, tp) are coupling 
and (total) potential functions of scalar fields, respectively. Also, A((f>, tp) and B((j), ip) are two typical 
functions coupled to the kinetic terms. We assume homogeneous scalar fields, that is 4> = <p(t) and 
V> = Y>(£), i n the following FRW metric 

ds ? = _ N 2 {t)dt 2 + a 2 {t) l^_ + r 2 d tf\ (2) 

where N(t) is a lapse function, a{t) is a scale factor and k specifies spatial geometry of the universe. 

We restrict our considerations to a non-interacting conformally scalar field model. This ansatz 
is obviously restrictive. However, a general reason for selecting such scalar fields is for simplicity. 
Indeed, it allows exact solutions in simple cases, as those will be discussed in this work and are rich 
enough to be useful as a probe for significant modifications that NCG introduces in the classical and 
quantum cosmologies. Thus in this case, one can set F(cj), ip) = /(</>) + f(ip) where / is a function of 
its argument as fix) = V(^ K ) — £x 2 /2- Also, k = SttG/c 4 and £ is a non-minimal coupling parameter 
with an arbitrary value that represents a direct coupling between the curvature and the scalar fields. 
The case £ = obviously is the minimally coupling situation; however in this work, we consider the 
conformal coupling case, i.e. £ = 1/6, and employ the units % = 1 = c and k = 3 (i.e. G = 3/8ir). 
We also, in general, consider the scalar fields that are not subject to any potential^ and assume more 
simple cases with A((j),ip) = 1 = B((p,ip). However, at the end of this section, we will investigate the 
particular case of non-zero constant values of the potential function (i.e. the cosmological constant). 

Now, based on these assumptions, substituting metric (2) into action (1) yields the Lagrangian 
density 

C = ( kNa -f](l- *±£) + (* + **)% + ^(i* + i?) , (3) 



1 A vanishing (total) potential, V(4>,'i]S), can also be achieved by non-interacting (inner) potentials, e.g. V{(j>,'ip) — 
U{4>) + W(ip) = A + (—A) = 0, where A is the cosmological constant. Such a case is an important one in models with 
double scalar fields [9]. 



where a total time derivative term has been neglected. By rescaling the scalar fields as 

x = a<j)/v2 and y = atp / v2 , (4) 

Lagrangian (3) reads 

C = kNa- a -£ + a{±2 + y2) - kN{x2 a +y2) . (5) 

Then, the corresponding Hamiltonian is 



n = N 



-,2 _i_ J2 



Pt , Px+Py f „ , /c(x 2 + y 2 ) 



4a 4a 



A;a + 



(6) 



where p a , p x and p y are the corresponding canonical conjugate momenta. For the conformal time 
gauge selection, namely N = a, one gets 



n = _Pl + £+A _ ka 2 + k(x 2 + y 2y 



(7) 



Therefore, the Hamilton equations are 



a = {a,l-L} 



1 

Pa = {Pa,%} = 2fea, 

x = {x,%} = -p x , 
Pi = {Px,"H} = -2fex. 

y = {y,K} = -p y , 

P y = {Py,!!} = -Iky. 



(8) 



Solutions of equations (8), with the Hamiltonian constraint % ~ 0, depend on the curvature index; 
actually their solutions are 

a(t) = A\ cost + A2 sini, x(t) = B\ cost + B2 sin i and y(t) = Ci cost + C2 sin t , 
/- - 1 : I (9) 

with constraint : 2{A\ + A%) = Y,l=i(Bf + Cf) , 

a(t) = A^e 1 + A 4 e~ t , x(t) = B 3 e f + B 4 e _i and y(t) = C 3 e* + Qe - * , 
A' --!:<! (10) 

with constraint: 2A3A4 = -B3-B4 + C3C4 , 



A; = 



a(i) = A 5 t + A 6 , x(t) = B 5 t + B 6 and j/(t) = C 5 t + C 6 , 



with constraint : 2A§ = B\ + Cf , 



mi 



where Aj's, Sj's and Cj's are constants of integration. As it is obvious, if one of the scalar fields 
vanishes and/or be equal to each other, i.e. <fi = ^(= x/Vty, all equations will lead to those derived 
in the case of one scalar field in our previous work [24] . 

We will compare these solutions with their NC analogs in the next section. 



2.2 The Noncommutative Phase Space 

In the classical physics, Noncommutativity is achieved by replacing ordinary product with the 
so-called Moyal product, shown by the * notation. This associative product on 2Z-dimensional phase 
space is defined between two arbitrary functions of general phase space variables, namely ( a = (<? i ,p ? ) 
for i,j = 1, • • • ,1, as [15] 



(f*9)(Q =exp 
such that the symplectic structure is 



2 a b 



/(C1MC2 



(12) 

Ci=C2=C 



<«*>-(V-, % ')' <13) 

where a,b = 1, 2, • • • , 11. The matrix elements a a b are assumed to be real, Oij and (3ij are antisym- 
metric, and Gij (which can be written as a combination of Oij and fy) is symmetric. The modified 
(or the a-star deformed) Poisson brackets are defined, in terms of the Moyal product, as 

{f,g}a = f*g-g*f- (14) 

Hence, the modified Poisson brackets of the phase space variables are 

Uii Qj}a = Oij , {li,Pj}a = $ij + <?ij and {pi,Pj} a = Pij ■ (15) 

For classical evolution, one starts by proposing that the NC nature of the minisuperspace variables 
can be encoded in a deformation of the Poisson structure. The main advantage of this approach is 
that the Hamiltonian does not need any modification. Hence, for our model, we require that the 
algebra of the minisuperspace variables obeys relations (15) in terms of the Moyal product defined 
in (12). However, to introduce such deformations, it is more convenient to consider a linear and 
non-canonical transformation on the classical phase space as (see e.g. Refs. [25, 26]) 

1 1 

x'i = Xi - -OijP 3 and p' i =p i + -f3ijX 3 . (16) 

Then, if the variables of classical phase space obey the usual Poisson brackets, i.e. {xi,Xj} = = 
{Pi,Pj} and {xi,pj} = Sij, the primed variables will yield 

{x'i, x'j} = Oij , {x'ijp'j} = S^ + aij and {p'^p'j} = fy , (17) 

with aij = —9yi(3j\[5 kl /4:. Consequently, to consider the noncommutativity effects, one can work 
with the ordinary Poisson brackets (17) instead of the a-star deformed Poisson brackets (15). In- 
deed, transformation (16) allows an extension of the commutative classical mechanics to the NC 
one. In the geometrical language, the usual Poisson brackets are mapped onto the modified Poisson 
brackets through transformation (16) and vice versa. It is evident that, for a compatible extension, 
transformation (16) must be invertible, and this imposes a condition on the NC parameters which we 
will specify for our model in below. 

Furthermore, let % = ~H(xi,Pi) be the Hamiltonian of a system including the commutative vari- 
ables; we shift the canonical variables through (16) and assume that the functional form of the 
Hamiltonian in the NC case is still the same as the commutative one, i.e. 

"H nc = "H(x-,p'J =n(xi- -O i jp j ,p i + -(3ijX j J . (18) 

This function is also defined on the commutative space and, obviously, equations of motion for 
unprimed variables are x % = &H nc /dpi and p 1 = —&H nc /dxi. Evidently, the effects due to the 
noncommutativity arise by terms including the parameters Oij and fyj. 
Now, in our model, the following notations are adopted 



(x\x 2 ,x 3 ) = (a,x,y) and (p 1 ,^ 2 ^ 3 ) = {p a ,Px,P y )- 

And, in this work, we assume that the only non-zero NC parameters are 2 

9 23 = 9 > and (3 23 = 4(3 > ; (19) 

hence <7 22 = a 33 = 9/3, where 9 and (3 are constants. Actually, in accordance with (16), we consider 
the following transformation 

a n 

x -> x - -p y , p x ~>Px + 2fiy , y ->■ y + 2^ and Py^Py- 2 P X - ( 20 ) 

Transformation (20) dictates that it is inverted when its determinant is not zero, that is when 0/3 ^ 1. 
Then, by making transformation (20) in Hamiltonian (7), the NC Hamiltonian is 

^nc = -^ + ^^(^+^)-fca 2 + (fc + /3 2 )(x 2 + y 2 ) + (^ + /3)( OT;c - a ;p J/ ), (21) 

and equations of motion are 

a = {a, U nc } = --zPa » Pa = {Pa, "Hnc} = 2fca , (22) 

and 

x = {x, n nc } = ^(i + ke 2 ) Px + (ke + p) y , 

Px = {Px,Hnc} = -2{k + (3 2 )x + {k9 + (3)p y , 

V = {V, ^nc} = \( l + k ° 2 )Py - (kd + (3)x , 

Py = {Py, «nc} = -2(fc + /3 2 )y - (fe0 + /3) Pa; . (23) 

Equations (23) show that, in general, the motion equations of the scalar fields are coupled in the NC 
case, and these equations reduce to the commutative equations (8) for e = = /3, as expected. Also, 
in the purposed model, the noncommutativity does not affect the time dependence of the scale factor 
and its solution is the same as the commutative case. 

The motion equations of the scalar fields, after eliminating momenta variables in (23), are 

x = -k(l-ei3) 2 x + 2(ke + l3)y and y = -jfe(l - 60) 2 y - 2{k0 + p) x, (24) 

with the Hamiltonian constraint 

x 2 + y 2 = constant if 1 + k6 2 = 3 

or 

(x 2 + y 2 ) + /c(l-6>/3) 2 (x 2 + y 2 ) = constant if 1 + k6 2 / 0. 

Space-time geometries with k = 0, 1 only satisfy the latter constraint for any 9, while k = — 1 geometry 
fulfills the former constraint with 9=1 (and consequently, because of the inversion condition, /3 7^ 1), 
and the latter one with 6^1. In general, solutions of (24) depend on the sign of a quantity defined 
as A = (l + k9 2 )((3 2 + k). 

For k = 0, 1, the sign of A is always positive, and thus, the real solution of equations (24) can be 
written as 

k = 0,l:{ x( ; t) = A A COSUJlt + ^ COSUJ2t (25) 



2 Note that, noncommutativity between x and y while they commute with the scale factor is completely consistent 
with noncommutativity between the original scalar fields, <j> and i\), while they also commute with the scale factor. 
3 The condition f + k6' 2 = is possible only when k — — 1, and hence 9 — 1. 



where A and B are constants of integration subject to the corresponding Hamiltonian constraint, and 
Ul = y/A-kd-P and uj 2 = -y/K - kO - (3. (26) 

Obviously, the corresponding terms in (25) are in tt/2 phase difference. 

In the case of k = 0, the NC parameter 9 does not actually appear in equations (24), and hence, 
in a spatially flat FRW universe, the scalar fields motion equations are affected only by the NC 
parameter (3. Besides, one gets uj\ = and oj 2 = —2/3, and thus solutions are similar to the motion 
of a free harmonic oscillator. However, the time dependence of solutions are modified from linear in 
the commutative case, to periodic in the NC case with the period of ir//3 which allows the model to 
be adjusted more easily with the observational data. 

For k = 1 geometry, the commutative solutions (9) are simple harmonics with the period of 2ir, 
whereas, the NC solutions (25) in general are not, though they still oscillate between two limits. 
Incidentally, solutions (25) are symmetric with respect to the NC parameters, and results do not 
vary by replacement of their roles. Besides, as ui\ ^ ui 2 , each solution of (25) can be considered 
as a forced (driven) harmonic oscillator, such that if the ratio of cj 2 /oji is a rational fraction, then 
solutions will be periodic (as in the Lissajous figures) with angular frequency given by the greatest 
common divisor of uj\ and uj 2 . Otherwise, solutions are non-periodic and never repeat themselves. 
A general condition that picks periodic solutions has been obtained in the appendix (supplemented 
with a diagram). Meanwhile, in the following special example, we provide a better insight about this 
situation. 

Without loss of generality, let us consider the special case ,5 = 0, with A = —B = D/2, for which 
solutions (25) read 



x(t) = D sin (9t) sin (y/l + 6 2 t) 
y(t) = D cos (9t) sin (y/l + 9 2 t) , 



with initial conditions x(0) = = y(0), x(0) = and y(0) = D\Jl + 9 2 . The condition for the 
periodic solutions is 



where m and n 2 are positive integers and obviously n 2 > n\. By solving this equation in terms of 9, 
one gets 



\ (n 2 - rai) /^n x n 2 = \ \y/n 2 jn\ - y/ni/n 2 ) ■ 



This relation provides different values of 9 parameter in terms of two integers, for which solutions 
(27) can be periodic. It is constructive to plot 9 in terms of r = n 2 jn\ as a continuous quantityf and 
this has been illustrated in Fig. 1 for range I < r < 80. As it is obvious, d9/dr is always positive, thus 
the function 9{r) changes monotonically with r. Also, when r has very large values, the rate d9/dr 
approximately decreases as 1/y/r. For the smallest allowed values of integers, i.e. n 2 = 2 and n\ = 1, 
we get the smallest value of 9 parameter, i.e. 9mm = 1/2 V2 for the periodic solutions. In solutions 
(27), the terms Dsin(9t) and D cos(9t) are envelopes of the corresponding curves. These types of 
oscillations, due to the periodic variation of the amplitude, are usually called beats. The well-known 
example of such oscillations is the simple harmonic one, where a driven force causes mechanical beats. 
Also, this situation can be simulated in the acoustical systems that produce a sound effect, known as 
beating [27], when \ui 2 \ — \uj\\ = 29 is sufficiently small and the terms including 9 2 in (27) are ignored. 
As a sample illustration of the beating effect for these cosmical oscillators, we have graphed a plot of 
x(t), equations (27), with its envelopes ±Dsm(0t), in Fig. 2 (right) for numerical values 9 = 0.1 and 
D = 1. As the figure illustrates, beats are more obvious when the envelopes are drawn. 



4 Note that, 712/ni is a rational number greater than unity, indeed, r = 1 yields 9 — which resumes the commutative 
case. 
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Fig. 1: The NC parameter 9, in k = 1 case with /3 = 0, as a function of r = 112/711. The points whose r are rational 
fractions correspond to the periodic scalar fields, and other points correspond to non-periodic ones. 





Fig. 2: The NC case for k = 1: periodic field (left) with 9 — 2 = /3, beating effect (right) with 
including the cosmical oscillator (solid line) and envelops (dashed lines). 



0.1 and p = 



y 2 




Fig. 3: Non-periodic field in the NC case for k — 1: with 8 = 2 and f3 — 0.1 (solid line), 
and with 9 = 0.2 and /3 - 0.01 (dashed line). 



When the ratio of 0J2 and u>\ is not a rational fraction, solutions are non-periodic but, as periodic 
cases, their behaviors still depend on the values of the NC parameters. For example, if one constructs 
a plot with numerous or a few relative extremum in a given time interval, then iterative drawings 
will indicate that the separation between the high and low points increases when the NC parameters 
tend to smaller values. This property is intensified for values less than unity, which is illustrated in 
Fig. 3 for y(t) in equations (25) with constants A = 2 and B = 1. From this point of view, the NC 
solutions have particular preference with respect to the corresponding commutative ones. 

For k = —1 geometries, one has A = (1 — 9 2 ){f3 2 — 1) that can be either positive, or negative 
or zero depending on different choices of the NC parameters. That is, when (9 > 1 and (3 < 1) 
or (9 < 1 and (3 > 1), A is positive and hence one gets solutions of type (25). However, in this 
case, the NC parameters have upper and lower bounds, and there is more restriction on finding the 
periodic solutions than in k = 1 geometry. On the other hand, when A is negative and 6^/3, the 
real solutions of equations (24) are 



x(t) = (C sinh y^At + D cosh V^At) cos(6 - (3)t 
y (t) = (C sinh y^At + D cosh V^At) sin(0 - /3)t, 



(28) 



where C and D are constants of integration subject to the corresponding Hamiltonian constraint. 
Solutions (28) are constrained to the condition that both NC parameters simultaneously have a lower 
bound, namely 9 and j3 > 1, or an upper bound, 6 and f3 < 1. Also, in the case k = 1, oscillations 
of the scalar fields in (28) have the phase difference of tt/2. These solutions, in contrast to their 
commutative analogs, are oscillating with a time hyperbolic amplitude that depends on the NC 
parameters and gives enough room for better adjustments. For instance, one can change the time 
interval between two successive zero points of the scalar fields, for the interval is tt/\6 — (3\. Besides, 
increasing and decreasing amplitudes in (28) with respect to the time depend on initial conditions. 
As an example, if one chooses the initial condition C = —D, then solutions of (28) will oscillate with 
decreasing amplitudes. Such solutions are similar to the damped harmonic oscillators with amplitudes 
proportional to exp(— \/— At) as envelopes. The main characters of such an oscillator, namely the 
time decay and the natural frequency, can be described in terms of the NC parameters as r = 1/y— A 
and ujq = ^oj\UJ2 = |1 — 0(3\, respectively. Another interesting point is that, if one assumes 6 and 
(3 > 1 then, solutions will damp quickly. Inversely, to possess more lately damped oscillations, the 
best choice of the NC parameters in the range 9 and (3 < 1 is when one assumes f? — ^ 1 — and (3 — > 0, 




Fig. 4: Damping scalar field in the NC case k = — 1: late damping with 9 = 0.999 and /3 = 0.001 (solid line), and 

quick damping with 9 — 1.5 and /3 = 1.05 (dashed line). 



in which each of the scalar fields has a maximum number of oscillations before the complete damping 
occurs. The diagram of such an oscillation is plotted in Fig. 4 for the x(t) of solutions (28) with D = 1 
for numerical values (6 = 1.5 and j3 = 1.05) as quick damping, and for (8 = 0.999 and j3 = 0.001) as 
late damping. 

When 8 = j3 va. k = — 1 geometry, solutions of (24) are decoupled hyperbolic time functions with 
the coefficient (1 — 8 2 ) in the exponent, which is the only difference between the NC and commutative 
solutions. The case A = is also possible when 8 = 1 or when /3 = 1, where solutions (28) are again 
periodic with the period of 2tt/\/3 — 1| or 2ir/\6 — 1|, respectively. Note that, the case k = — 1 with 
8 = 1 resembles the constraint condition 1 + k8 2 = with arbitrary f3, and however for f3 = 1 one 
gets trivial constant solutions, but this value is not allowed. 

Non-Zero Constant Potential 

Let us obtain the equations of motion for the scale factor when one considers a non-zero constant 
value for the potential function (i.e. a cosmological constant) in action (1). Thus, suppose the 
potential function is V((p, -0) = A. In a manner similar to the free-potential case, one gets the 
Hamiltonian 



n _ Pa , Px+Py 

4 4 



ka 2 + k{x 2 + y 2 ) + Aa 4 



(29) 



in the commutative case, and 



W. 



- P -f + —-—(p 2 x +p 2 y )-ka 2 + (k + p 2 ){x 2 +y 2 ) + {k8 + P)(yp x -xp y ) + ka\ (30) 



in the NC one. It is obvious that the Hamiltonian equations for the scalar fields, derived from (29) and 
(30), are the same as those obtained in the corresponding free-potential case, namely equations (8) 
and (23), respectively. But, the Hamiltonian equations of the scale factor, for both the commutative 
and NC cases, are modified in the same manner as 



a = — - 



Pa 



and 



p a = 2ka - 4Aa 3 . 



(31) 



Eliminating the momentum variable gives 



a = 2AcT - ka. 



(32) 



m 



By integrating, one gets 

a 2 = Aa 4 - ka 2 + t> 2 , (33) 

where Vq is an integration constant. Solution of equation (33) can be written in terms of the Jacobi 
elliptic functions, that is 

a(t) = — sn (at, m) , (34) 

a 



where the initial condition a(0) = 0, 



-1/2 



a 2 = - I k + \jk 2 - 4 Av 2 J and m = [(a 2 /Av 2 ) - lj . (35) 

The Jacobi function, on the right hand side (34), behaves as a sine function when a is a real and 
m is either a real or a pure imaginary number? Therefore, for negative potentials in an arbitrary 
geometry, and also for positive potentials in the k = 1 geometry when 4Av 2 < 1, the scale factor 
behaves periodically as a sinusoidal function. On the other hand, for positive potentials in k = 0, — 1 
geometries, and in the k = 1 geometry when 4 Av 2 > 1, the Jacobi function has a complex value, and 
hence these situations are not allowed. 

In the limit A — > for k = 1, one gets m — > 0, a 2 — > 1 and sn(t,0) behaves as sin(t). Hence, we 
get back to the free-potential solution (9), as expected. However, in such a limit, for k = —1,0, as 
a — > 0, solution (34) is not valid, and hence, one must consider equation (33). Now, by regarding 
the limit A — > in equation (33) for k = —1,0, and solving the resulted equation, one again gets the 
free-potential solutions (10) and (11), respectively. 

3 The Quantum Model 

Although the effects of NC cosmology are mostly desired and important in classical approaches, 
but still it is instructive to investigate its quantum counterparts. In particular, if a universe has been 
commenced by a big bang or from a very very tiny scale, then it would be plausible that quantum 
behaviors should have significant implications in its evolution and cosmology. Thus, in this section 
we proceed to quantize the cosmological model given by action (1) for free-potential, such that 
the canonical quantization of the phase space leads to the WD equation, Ti^ = 0, where % is the 
Hamiltonian operator and ^ is a wave function of the universe [29] . We employ the usual canonical 
transition from classical to quantum mechanics via the generalized Dirac quantization of the Poisson 
brackets to quantum commutators, i.e. {} — > —i[]. Then, as the classical approach, we investigate 
the general properties of the wave function in the commutative and NC frames of the quantum model 
in the following subsections. 

3.1 The Commutative Quantum Cosmology 

As usual, the operator form of Hamiltonian (7) can be acquired by the replacements p a —> —id a , 
p x — > — id x and p y —> —id y . Assuming a particular factor ordering, the corresponding WD equation 

is 



da 2 \ dx 2 dv 2 I V 



a 2 



*(a,x,y) = 0. (36) 



da 2 \ dx 2 dy 2 
In terms of the polar coordinates 

x = pcosip and y = psimp, (37) 



5 For properties of the elliptic functions, see, e.g., Refs. [28]. 
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equation (36) reads 



d 2 a -4ka 2 -(d 2 + ^ + %)+4kp 2 



V{a,p,<p) = 0. 



P P 

Let us choose a product ansatz as a solution of equation (38) , namely 

*(a,P,<P)=Ma)B(p)e 2 ^, 

where v is a real constant. By substituting (39) into equation (38), one gets 

A" + 4 (p - fca 2 ) A = 

and 

„2 D // 



(38) 



(39) 



(40) 



p 2 B" + pB' + 4 W 2 - fep 4 - v 2 )B = 



(41) 

where the prime denotes ordinary derivative with respect to the argument and p is a constant of 
separation. Solutions of equations (40) and (41) for real values of p, corresponding to the curvature 
index and boundary conditions, are 



k = : 



A-n (a) = C\ sin (2^~pa) + C 2 cos (2^/pa) 
B^ (p) = D x J 2v (2^pp) + D 2 Y 2v (2^pp) 



(42) 



and 



\(a) 



,-V2 



C z M_ iL _ 1 [2Vka 2 ) + C 4 W_ E _ 1 [2-Jka 2 

2vT'4 V / 2vT'4 V 



-1,1: < 



(43) 



B„ v (p) = p- 1 



D 3 M^_ „ (2^/kp 2 ) + D A W. 



2Vfc' 



2vTfc' 



2\fkp'< 



where J 2v and Y 2v are respectively the Bessel functions of the first and second kind, M Vt \ and W„\ 
are the Whittaker functions and, Cj's and D^s are constants. For positive (zero) curvature, the 
Whittaker (Bessel) term M V: \ (Y 2u ), in the classically forbidden region, is divergent? Thus, one can 
discard these terms and write the well-defined eigenfunctions of equation (38) as 



ipfMv(a,p) = A fl (a)B flv (p) oc < 



Ci$m.(2y/pZa) + C 2 cos (2^/pZ a) J 2l/ (2^/pp) 
{p^)~ l W E i (2a 2 ) We iV (2p 2 ) 
(p^r'M^ 1 (2ia 2 ) M_ t tL „ (2ip 2 ) 



for k = 0, 

for fc = 1, (44) 

for fc = — 1, 



where in k = — 1, for simplicity, we have written solution (44) only in terms of the Whittaker function 
Mjj^x- The wave packet corresponding to (44) is 



* 



/CO roo 
/ E^E^^e^^dpdv . 
-00 J — 00 



(45) 



where -E 1 ^ can be taken [16, 19] to be a shifted Gaussian weight function with constants b and c as in 
exp[— b(p — c) 2 ], and a similar expression for E v . 



^For the properties of the Whittaker and Bessel functions and their indices, see again Refs. [28]. 
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3.2 The Noncommutative Quantum Cosmology 

It is well-known [26] that, in the NC quantum mechanics, the original phase space and its sym- 
plectic structure are modified. That is, for the NC proposal of quantum cosmology, we assume that 
operators (variables) of the FRW minisuperspace obey a kind of deformed Heisenberg algebra like 
the ones in the NC quantum mechanics as [26] 



id 



K) ■ 



[oti,Pj] = i(Sij + Oij) and [pi,pj] = ifaj ■ 



(46) 



The notations and definitions are the same as in the NC classical model. This kind of extended 
noncommutativity maintains symmetry between the canonical operators, and yields the usual Heisen- 
berg algebra in the limit 0y and fy — > 0. As usual, this deformation can be redefined in terms of 
noncommutativity of minisuperspace functions with the Moyal product defined in (12). Thus, the 
corresponding noncommutative WD equation can be written by replacing the operator product, in 
equation (36), with the star product, namely "H * \P = 0. 

However, it is possible to reformulate equations in terms of the commutative operators with the 
ordinary product of functions if the new operators are introduced, again as our previous assumption 



yPy 



Px 



p x + 2(3y , y =y + -p x and p y = p y - 2(3x . 



(47) 



Clearly, if unprimed operators obey the usual Heisenberg commutators, then the non-zero primed 
operators will obey the deformed Heisenberg commutators (46) in the form 



[x',y']=i9, [x',fl x ]=i(l + 



[if, fly] and [fl x ,fl y ] = 4i(3. 



Therefore, operator transformation (47) can be regarded as a generalization of the usual quantum 
mechanics to the NC one. On the other hand, the inverse transformation of (47) are 



x = rj(x' + \fl) , p x = r](fl x - 2f3y') , y = rj(y' 



2Px) 



and p y = 7] (fly + 2/3 x') 



where r] = 1/(1 — 9/3). Consequently, one can go from the usual commutators to the deformed ones 
and vice versa provided that again 9(3 7^ 1. As a result, the original equation, employing the new 
operators, reads [30] 



1 



1 



H(xi,pi) * ^ = H[xi- -9ijfl,Pi + -fiijx 1 ] * 



0. (48) 

Hence, the noncommutative WD equation corresponding to the NC Hamiltonian (21) is 

~d 2 a - Aka 2 - (1 + k9 2 ){d 2 x + d 2 y ) + Ai{k9 + (3)(xd y - yd x ) + 4(fc + (3 2 )(x 2 + y 2 )] * nc (a, x,y)=0, 

(49) 
which in the polar coordinates (37) yields 



d n d[ 



2\„2 



d l a - Aka z - (1 + k9 z ) K + — + -f + &{k9 + /3)^ + 4(fc + /3> 



* nc (a, p, (p) = 0. (50) 



When 1 + k9 2 / 0, by using the ansatz (39), equation (50) is separable to 



and 



p 2 B" + pB' + 4 



A" + 4(/x - ka 2 )A = 0, 



p + 2v{k9 + p) 2 k + /3 2 



l + k9 2 



l + k9 2 



P 



B = 0, 



(51) 
(52) 



7 As will be seen, this condition should automatically be satisfied when the wave function is well-behaved. 
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where ^ is a constant of separation. The scale factor part of the wave function, equation (51), is 
similar to its commutative analog (40), as expected. The radial part of the wave function, equation 
(52), first reduces to its commutative analog, equation (41), when 9 = = (3, as again expected. 
Secondly, in the case k = with (3 = 0, equation (52) once again reduces to equation (41) even if the 
NC parameter 9 does exist. However, in k = 0, solutions of equation (52) do not depend on the NC 
parameter 9 at all. Namely in a flat FRW universe, the (3 parameter is the only responsible parameter 
for the NC effects. These properties are common with the classical model. 

Comparing equation (52) with the commutative analog (41) shows that the functional form of the 
radial part of the wave function (and hence, the whole wave function) is the same as the commutative 
ones provided that the coefficients p and k in equation (41) are replaced by 

M^ l + k02 and k ^ l + k9 2- W 

Therefore, for k = 0, 1 and k = —1 with 9 ^ 1, the NC eigenfunctions are in terms of the Whittaker 
functions similar to the commutative solutions (44). However, in the case k = — 1, there are bounds 
on the NC parameters. For a better illustration, let us write (53) as 

p ->■ o" = /J+ i_giT and k ->• K = -jz§z- 

Since this modified curvature index, K, can have positive, negative and zero values, the solutions 
again are similar to the commutative cases (42) and (43), where K = also requires (3 = 1. But, for 
K/0, the NC parameters must satisfy inequalities depending on the sign of K. Namely if K < 0, 
then there will be (8 and (3 < 1) or (9 and (3 > 1), and if K > 0, then (9 > 1 and (3 < 1) or (9 < 1 
and (3 > 1). For instance, by choosing inequality (9 and (3 < 1) when K < 0, the corresponding wave 
packet can be written in the form 

tf nc (a, p,<p) = / E li E v (py/a)- 1 M_ iE i hia 2 ) M°, v h^/K p 2 ) e 2iu *dpdu , (54) 

where the lower limits, in contrast to the commutative wave packet (45), are bounded by p + 2v((3 — 
9) >0. 

In the case 1 + k9 2 = 0, namely when k = — 1 and 9 = 1 (hence (3^1), with the separation of 
variables ^ nc = A(a)B(p, (p), equation (50) reduces to two differential equations. One equation is 
similar to equation (51), and the other one is 

i(l-(3)d v + (l-p 2 )p 2 ]B(p,<p) = C, (55) 

where C is a separation constant. As, in general, C has a non-zero value, equation (55) does not 
have a suitable solution for (3 = 1 (which itself is not allowed). For (3 ^ 1, one easily shows that its 
solution can be written as 

B(p, i p) = F(p)e^ 1+ ^^ + B /p 2 , 

where F is an arbitrary function and Bq = C j(l — (3 2 ). As this function is not well-behaved when p 
tends to zero, it is not allowed. 

4 The Noether Symmetries 

In this section, we employ the Noether theorem and explore the effects of noncommutativity on 
the underlying symmetries in the commutative frame. For this purpose, following Refs. [9, 31], one 
can define the Noether symmetry as a vector field, say X, on the tangent space of the phase space. 
In our model, it can be, in general, as 

X = A— B— C— — — — — — — (56) 

da dx dy dt da dt dx dt dy ' 



such that the Lie derivative of the Lagrangian with respect to this vector field vanishes, i.e. 

L X C = 0. (57) 

For simplicity, we assume that unknown functions A, B and C to be linear in terms of the configuration 
variables a, x and y. The time derivative d/dt represents the Lie derivative along the dynamical vector 
fields, which in our model is d/dt = ad /da + xd/dx + yd/dy. 

Now, in order to obtain constants of motion, let us rewrite equation (57) as 



dC dA dC 
da dt da 



dC dBdC 
dx dt dx 



+ 



dC dCdC 
dy dt dy 



By employing the Lagrange equation, dC/dq= dp q /dt, it reads 



dp a dA 

A— — p a 

dt dt F 



T ' B itt + it p > 



+ 



dC 



c dpy_ 
6 dt + dt py 



0. 



0, 



which yields 



— (A Pa + B Px + C Py ) 

dt 



Therefore, the constants of motion are 



Q = A Pa + Bp x + Cp 



;</• 



(58) 



(59) 



(60) 



(61) 



for different unknown functions of A, B and C. To obtain these functions, one can employ equation 
(58) or equation (59) which is more suitable in the Hamiltonian formalism. To manage this, one can 
write equation (59) in terms of the Poisson bracket, { , H} = d/dt, as 



+ 



A{p a , %} + B{ Px , %} + C{ Py , %} + 

Px + 



dA, „,_ dA ( nil dA. ni - 

— {a,H} + j-ix^} + j-{y,n} 
da dx dy 



Pa 



dB . ni . dB . nn dB . ni . 

-K-{a,n} + -7-{x,u} + -rr{y,n} 

da dx dy 



dC . nil dC . ni . dC . „,■ 
— {a,H} + -t-{x,U} + -T-{y,U} 
da dx dy 



p y = 0. (62) 



This equation, in general, gives quadratic polynomials in terms of the momenta with coefficients being 
partial derivatives of A, B and C with respect to the configuration variables. Hence, the expression 
identically is equal to zero if and only if these coefficients vanish, which lead to a system of partial 
differential equations for A, B and C. 

In the following subsections, we obtain such symmetries for the model in the commutative and 
NC cases. 



4.1 Symmetries in The Commutative Frame 

In this case the Hamiltonian is given by relation (7); hence by substituting the corresponding 
Poisson brackets into equation (62), one gets 



2k(aA -xB- yC) 



+ - 



dA 

da, 



Pa + 



dB 



dC 



dx Px+ dy P ' y 



1 

+ 2 



dA dB\ 1 /dA dC 

dx da J 2 V dy da 



1 (dB dC\ 

P*Py +ohr + iu VxVv 



o. 



2 V dy dx J 
Let us first treat spatially non-flat geometry, k ^ 0, for which equation (63) leads to 

0, 
0. 



aA-xB - yC 
dA _dB _ dC 
da dx dy 

dA dB\ (dA dC\ (dB dC 

+ 



dx da 



dy da 



dy dx 



0. 



(63) 



(64) 
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A general solution of constraints (64) is 

A = c\x + C2y , B = c\a — c 3 y and C = 02a + c 3 x , (65) 

where Cj's are three constants of integration. For an easier representation of the commutation relations 
between symmetric vectors, let us change x — > —x and y — > —y in solution (65). Thus by (61), three 
independent constants of motion are 

Q\ = aPx ~ xp a , Q2 = apy - yp a and Q 3 = yp x - xp y , (66) 

which are the well-known angular momenta about the configuration variables. The corresponding 
symmetric vectors are 

d d . d . d 

Xi = a- x— + a— - x— , 

ox oa ox oa 

d d d d 

oy oa oy oa 

d d . d . d 

X 3 = y^. x— + y— - x— , 

ox oy ox oy 

(67) 

which satisfy the Lie algebra [A^X,-] = e^ ■ Xk , where e^ is the Levi-Civita tensor. 

For k = 0, it is clear from equation (63) that the symmetries can be obtained from the last two 
constraints of (64). Thus, the corresponding solution is 

A = d\x + d,2y + d 3 , B = d\a + d^y + d§ and C = o^a — d^x + d$ , (68) 

where d^s are six constants of integration. Therefore, the six independent constants of motion are 

Qi=Pa, Q2 = Px , Q3 = Py , 

(69) 
Qa = VPx - xp y , Q 5 = ypa + apy and Q e = xp a + ap x . 

Incidentally, in this case, a glance at Hamiltonian (7) shows that all configuration variables are cyclic 
and consequently their corresponding momenta are constants of motion, i.e. Q\ to Q3. In addition, 
Hamiltonian (7) for k = (actually the Lagrangian x 2 + y 2 — a 2 ) is invariant under rotation about 
the a-axis. Thus, the angular momentum about this axis is conserved, i.e. Q^. 
The corresponding symmetric vectors, in the flat geometry, are 

x - d x - d x - d 
oa ox oy 

d d d d 

ox oy ox oy 

d d d d 

Xs = yiT + a 7T + y~pF + a !T ' 
oa oy oa oy 

d d d d 

Xq = x— + a— + x— + a— , 
oa ox oa ox 



which satisfy 
















[Xi,X 2 ] 


= 0, 


[Ai,A 3 ] 


= 0, 


[A!,A 4 ] 


= 0, 




[Xi,X 5 ] 


= x 3 , 


[Ai,A 6 ] 


= x 2 , 


[A 2 ,A 3 ] 


= 0, 




[x 2 ,x 4 ] 


- -A3, 


[x 2 ,x 5 ] 


= 0, 


[A 2 ,A 6 ] 


= x u 




{x 3 ,x 4 } 


= x 2 , 


[x 3 ,x 5 ] 


= Xi, 


[A 3 ,A 6 ] 


= 0, 




[x 4 ,x 5 ] 


= —X&, 


[A4, A 6 ] 


= x 5 , 


[X 5 , A 6 ] 


= A 4 . 



(70) 
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4.2 Symmetries in The Noncommutative Frame 

Now, let us find out which of the above symmetries survive in the NC case. Here, the Hamiltonian 
is given by relation (21), and the required Poisson brackets are given by equations (22) and (23). 
Substituting the corresponding Poisson brackets into equation (62) gives 



1 

+ 2 



2kaA - 2{k + (3 2 )(xB + yC) 
9i, I, n x (dB 2 dC 2 



+- 



[l + w 



dA 



dx 
+{k9 + /3) 



dB 
da 
dA 



PaPx + 



(i + ke 



dA 



dC_ 

dy da 
dB dB 



PaPy + -(l + k9 )^— + — ]l>vPll 



dA\ ( dB dB \ ( dC dC 

dx dy J \ dx dy J \ dx dy 



0. (71) 



Obviously, equation (71) yields extra restrictions on A, B and C with respect to the commutative 
case (63). 

When k9 + j3 ^ Of by putting each coefficient in equation (71) equal to zero, one gets 

kaA - (k + (3 2 )(xB + yC) = 0, 
dA _ dB_ _ d£ _ 
da dx dy 



dA 



dB 
dx da 
dA dB 



M 7n2 ,dA dB ,., , n2 ,dA dC dB dC 
'! + M 2 )-5- - -£- = (1 + kO 2 )— - — = — + — 



y- 



dx dy 

Hence, one obtains the solution 

and 



y 



dx 



dy da dy dx 

dB dc dc n 

'ir- c = yjr- x jr + B 

dy dx dy 



0. 
■ 0. 



B = h y, 



A 



for 

for 



C 



k = 
k^O. 



-hnX 



(72) 



(73) 



(74) 



where ho and Aq are integration constants. Thus, when k ^ 0, the only constant of motion is 

Qnc = VPx ~ xp y , (75) 

with the symmetric vector 



X,, 







d 



d 



d 



dx dy dx dy 



(76) 



which is an especial case of commutative solution (65) with initial conditions c\ = = c 2 - For k = 
(with /3 ^ 0), in addition to Q nc , we have an another constant of motion, namely Q a = p a , which 
corresponds to the symmetric vector X a = d/da with the trivial Lie algebra LY nc ,X a ] = 0. This case 
also restricts commutative solution (68) with initial conditions d\ = d 2 = d§ = d§ = 0. Therefore, in 
the k6 + (3 t^ case, the noncommutativity reduces the number of symmetries to one for k ^ and 
two for k = 0. 

When k9 -\- (3 = 0, the last row of constraints in (72) is omitted. In this case, one choice is k = — 1 
and 9 = (3, where a non-trivial (i.e. 9 = j3 ^ 0) linear solution exists if and only if 9 = \[2 = f3. 
Hence, by the sign change of functions, i.e. B — > —B and C —> —C (or equivalently A — > —A), the 
solution is the same as the commutative solution (65), which again yields the constants of motion and 



The equality k9 + j3 = is possible when (k — — 1 and 9 = j3) and or (k = = /3 and any value of 9). 
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symmetric vectors (66) and (67). However for k = — 1 and a special value of 9 = 1 = j3, it gives one 
solution (i.e. (75)), but this case is not allowed due to the inversion condition 9/3 ^ 1. For another 
choice k = = f3 with any value of 9, it gives the commutative solution (68). That is, even though the 
noncommutativity is still present, the number of symmetries and constants of motion do not change 
with respect to the corresponding commutative case. 

Note that, in all cases, irrespective of whether the noncommutativity exists or not, and for any 
value of the curvature index, the angular momentum about the a-axis is conserved, as expected. 
Besides, in the above considerations, we have not, in general, specified numerical values of the NC 
parameters. 

5 Conclusions 

We have carried out an investigation for the role of NCG in cosmological scenarios, based on a four- 
dimensional free-potential (multi) scalar-tensor action of gravity, by introducing a NC deformation 
in the minisuperspace variables. The phase space is generated by two non-interacting conformal 
scalar fields plus the scale factor with their canonical conjugate momenta. The scalar fields are non- 
minimally coupled to geometry whose background is the FRW metric, where the conformal time gauge 
evolutions have been studied. The noncommutativity has been introduced only between the scalar 
fields and between their canonical conjugate momenta via two NC parameters 9 and /3, respectively. 
The investigation has been carried out for this toy model by means of a comparative mathematical 
analysis of the time evolution of the dynamical variables in the classical level and of the wave function 
of the universe in the quantum perspective, both in the commutative and NC frames. We have paid 
more attention to the outcome of results and have looked for the relations, including ranges and 
values, among the NC parameters for which particular or allowed solutions exist. 

In general, we have shown by this toy model that the purposed noncommutativity has important 
implications in the evolution of the universe, however, does not affect the time dependence of the 
scale factor, i.e. its solution is the same as the commutative case, as expected. Also, we have found 
that one of the particularity of the NC parameter in the momenta sector, i.e. /3, is in the spatially 
flat FRW universe, where it is the only responsible parameter for the NC effects in the classical and 
quantum frames. 

In the classical model, exact solutions have been obtained. One of the important aspects of the 
NC solutions is that they can be regulated with both NC parameters. For example, these parameters 
can be employed to adjust the time dependence of solutions with the experimental or observational 
data. A distinguished feature of the noncommutativity effects, which we call a cosmical oscillator, is 
that the scalar fields behave similar to (or can be simulated with) the three most important harmonic 
oscillators depending on three values of the spatial curvature. These are the free, forced and damped 
harmonic oscillators corresponding to the flat, closed and open universes, respectively. In the flat 
universes, the time dependence of solutions are modified from linear in the commutative case to 
periodic in the NC frame. In the closed universes, if the ration of frequencies of the scalar fields is a 
rational fraction, then solutions will be periodic. This condition restricts the NC parameters. When 
this ratio is not a rational fraction, the solutions are non-periodic but their behaviors still depend 
on the values of the NC parameters. A plot with numerous (or a few) relative extremum in a given 
time interval indicates that the separation between the high and low points increases when the NC 
parameters tend to smaller values, and this property is intensified for values less than unity. From 
this point of view, the NC solutions have particular preference with respect to the corresponding 
commutative ones. Furthermore, in the k = 1 case, the solutions are symmetric with respect to 
the NC parameters, and the results do not vary by replacement of their roles. Also, when the NC 
parameters are small, the cosmical oscillators have analogous effects with the familiar beating effects 
in the sound phenomena. For a better view on this situation, an example has been illustrated in the 
text. In the open universes, there are upper and lower bounds on the NC parameters. The quick 
and late damping of this case can be adjusted by the NC parameters. For example, for 9 — > 1~ and 
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(3 — > 0, the scalar fields have maximum number of oscillations before the complete damping occurs. 

We have also shown that the existence of a non-zero constant value of the potential function (i.e. 
the cosmological constant) does not change the time evolutions of the scalar fields, but it modifies 
the time dependence of the scale factor in a same manner for both the commutative and NC frames. 
Indeed, we have obtained that in all allowed conditions the scale factor behaves as a sinusoidal 
function. 

In the quantum model, the exact solutions for the wave functions of the universe have also been 
obtained. The scale factor part of the wave function is similar to its commutative analog, as expected. 
One still expects that when the noncommutativity effects are turned on in the quantum scenario, 
they should introduce significant modifications in the scalar fields. However, an interesting feature 
of the well-behaved solutions is that the functional form of the radial part of the NC wave function 
is the same as the commutative ones within the given replacements of constants caused by the NC 
parameters, although these replacements in turn may cause drastic effects. For example, the curvature 
index is modified, and in open universes, the allowed NC wave functions impose bounds on the NC 
parameters. 

Finally, we have employed the Noether theorem and have explored the effects of noncommutativity 
on the underlying symmetries in the commutative frame. We have shown that for spatially flat 
universes, there are six Noether symmetries, and, in general, only two of them are retained in the 
NC case. In the special case k = = (3, all symmetries survive regardless of the existence of 9 
parameter. For non-flat universes, there are three Noether symmetries in the commutative case, one 
of which is retained in the NC frame. However, in open universes, when the NC parameters have the 
values 9 = y/2 = (3, the number and general form of symmetries do not change with respect to the 
commutative frame. The only difference is related to the sign of symmetries. Conservation of the 
angular momentum about the scale factor axis is a common face between the commutative and NC 
cases, as expected in the purposed noncommutativity. 

We should emphasis that the scalar fields solutions are given after rescalings of the original fields, 
relations (4), and in the conformal time. Consequently, the functional forms of the solutions are not 
as simple when translated to the original fields. However, any detection should be performed in a 
reference frame, and there are debates on the physical frame in the cosmological contexts. Besides, the 
model is to be taken as a toy model only that still provides a valuable contribution to the assessment 
of the implications of the NC geometrical deformations of the phase space upon the dynamics of the 
cosmological model envisaged. However, a more realistic NC cosmological model may be achieved 
when one also involves the noncommutativity of the scalar fields with the scale factor, where the value 
of its time derivative in the form of the Hubble parameter can be determined from observations. 

Appendix 

We first obtain a general relation among the NC parameters when the solutions (25) are periodic. 
Then, we treat a specific periodic solution as an example. 

As mentioned in the text when k = 1, solutions (25) play as the forced oscillators, and the 
existence of periodic solutions is provided when the ratio u^A^i is a rational fraction. In this case, 
the A parameter, by its definition, is A = (9+(3) 2 + (9(3— l) 2 ; hence by the transformation reversibility 
constraint 6(3 ^ 1, one always has vA > 9 + (3. Therefore, by definitions (26), we can take u)\ = cn\ 
and UJ2 = —cri2, where n\ and n 2 are two positive integers with n 2 > m, and c is an arbitrary 
non-zero positive constant. Hence, again, definitions (26) give 

9 + l3 = -(n 2 -ni) and 9(3 = 1 ± c^/ralra^ > . (A.l) 

These relations impose a firm restriction on the values of the NC parameters (or equivalently on 
allowed positive integers). Indeed, the condition for periodic scalar fields is that, the values of the 
NC parameters must be such that there exist two positive integers satisfying relations (A.l). As 9 

1 o 



and (3 are real roots of the quadratic equation X 2 — (9 + (3)X + 9(3 = 0, relations (A.l) lead to 

(n 2 -n!) 2 > (16/c 2 )(l±c^T^)- (A.2) 

The above result implies that if one finds a pair of integers satisfying inequality (A.2), then the values 
of the NC parameters will be calculated from equations (A.l). 

Let us treat a specific example for the above considerations. Consider the equality in relation 
(A.2), which yields 9 = f3, which also by the constraint 9(3 ^ 1 gives 9 ^ 1 (note that we have 
assumed 9 > 0). Also, by definitions (26), we have u>i = (1 — 9) 2 and uji = —(1 + 9) 2 . Thus, by 
imposing the periodic condition as 

= I ) = rational fraction = m > 1, (A. 3) 

wi \l-9J 

it leads to 

9 = > 1 or 9 = < 1 , (A.4) 

m — 1 m + 1 

where m is a real number greater than one? In the limit m — > 1, one solution of (A.4) yields # — )■ 0, 
as expected; however, the other one gives infinity that is not physically accepted. Also, in the limit 
m — > oo, we have 9 — > 1, which, in general, may be viewed as a transition from the forced harmonic 
oscillator in the k = 1 geometry to the simple one in the k = case. In another words, when the 
value of m slowly gets bigger and bigger, then the driven force is gradually removed away from the 
cosmical oscillators. The time evolution of scalar field x(t), equation (25), is depicted in Fig. 2 (left) 
for the numerical value m = 3 corresponding to 9 = 2, with superposition constants A = 1 and 
B = 1. The period of solutions is the least common multiple of distinct periods, namely 27r/(l — 9) 2 
and 2tt/(1 + 9) 2 . 
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